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Abstract 

Total computational costs of scientific simulations are 
analyzed between direct numerical simulations (DNS) and 
multiphysics simulations (MPS) for sound generation in 
musical instruments. In order to produce acoustic sound 
by a turbulent flow in a simple recorder-like instrument, 
compressible fluid dynamic calculations with a low Mach 
number are required around the edges and the resonator of 
the instrument in DNS, while incompressible fluid dynamic 
calculations coupled with dynamics of sound propagation 
based on the LighthilVs acoustic analogy are used in MPS. 
These strategies are evaluated not only from the viewpoint 
of computational performances but also from the theoretical 
points of view as tools for scientific simulations of compli- 
cated systems. 



1. Introduction 

Scientific simulations have been widely used to study 
complicated dynamic properties of physical and chemi- 
cal systems along with the recent development of high- 
performance computing environments. However, when we 
investigate a large, complicated system, which involves 
multiple scaled dynamics in space and time forming a dy- 
namically hierarchical structure in the phase space, the 
method of the direct numerical simulation (DNS) faces a 
serious obstacle. That is, relatively large resources are re- 
quired to achieve precise description of the whole behav- 
ior of the target system. This is just the case for which 
we may introduce an alternative approach based on multi- 
physics simulations (MPS). In MPS, 1) the target systems 
are divided into smaller parts described by relatively simple 
physics; 2) each separated component is simulated indepen- 



dently according to the individually governing equation; 3) 
communication among separated components including the 
dynamical interactions is achieved by using "coupled sim- 
ulation techniques", in which a mediator plays a role of an 
information center for MPS to control information among 
the extensive computational resources. 

Acoustic sound generation in musical instruments f4l is 
an interesting target of the numerical study, for which a 
careful consideration in modeling physical processes as 
well as program codes is required to to reproduce exper- 
imental results. Many works have been published using 
empirical instrumental models d [71 [51 and it has been re- 
vealed that complicated dynamics are induced by the cou- 
pling between airflows and generated sounds. Recent stud- 
ies on sound production in air-reed instruments, e.g. flute, 
recorder, organ pipe and so on, often simulate jet oscilla- 
tion by a resonator or turbulent dynamics with vortex shed- 
ding using high-performance computing methods |[4l[5l[T11l. 
However, those approaches seem to achieve no more than 
partial success. Thus, it is generally difficult to execute 
DNS calculations for sound generation since the calcula- 
tions of compressible fluid dynamics in three dimensional 
configurations are still too heavy in the usual computer en- 
vironments. One of the difficulties of the sound simula- 
tion can be found in difference of the energy scale between 
turbulent flows and radiated sounds, where the energy of 
the airflow is 10^ times larger than the sound generated un- 
der typical conditions in playing musical instruments. It is 
practically impossible to study such problems of multi-scale 
properties by using the DNS. Indeed, the larger resouce is 
required to calculate detailed behaviors in the parts of the 
smallest scales. MPS or multiscale simulations (MSS) seem 
to provide a breakthrough to solve those problems. Further- 
more, technical advantages of the MPS approach, which is 
also called hybrid methods, are generally known in compu- 
tational aeroacoustics |[T2ll . 



On the mechanism of the sound production by fast air- 
flow around objects, there have been many theoretical and 
numerical works on the interaction between fluid dynam- 
ics and acoustic sounds I T2l . The standard approach by the 
Lighthill's acoustic analogy f6l is widely used both in the 
recent reseaches based on the DNS strategy and the approx- 
imate studies using Large-Eddy Simulations (LES) of in- 
compressible fluid dynamics. 

In this paper, we apply the MPS technique to the sound 
production by turbulent flows in musical instruments, where 
a consistent coupling between the turbulent flow dynamics 
and the sound propagation is considered. This paper is or- 
ganized as follows. The basic theoretical and mathematical 
descriptions on the sound generation are briefly reviewed 
in section O Based on the theory, the coupled simulation 
system of the turbulent flow and the sound propagation is 
constructed in section O In section IH computational costs 
for DNS and MPS approaches are evaluated, and the con- 
cluding remark will be given in the last section. 
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Figure 1. Schematic picture of multiphysics 
simulation system for musical instruments 



2. Theories for sound generation 

The origin of sounds in flutelike instruments is usually 
considered as the jet oscillations, vortex shedding, and tur- 
bulences. In this section, we briefly review the basic theo- 
ries on these sound production in musical instruments. 

2.1. Lighthill's acoustic analogy 

The sound radiated from a finite region of turbulent 
flow is estimated by Lighthill's theory on the acoustic anal- 
ogy. This provides a general approach to compute the field 
where the acoustic pressure is generated by a function of the 
Lighthin tensorlOlllIl 

Tij = pUiUj + - Co'^p)Sij - CFij. (1) 

Under some assumptions, it is known that only the main 
term of T^j, i.e., pUiUj, is necessary. Then, the wave equa- 
tion of sounds for the density of the fluid is 

It is important that the Lighthill equation is an exact re- 
statement of the Navier-Stokes equations Ol. This is the best 
starting point for the multiphysics approach explained in the 
next section since the mathematically exact representation 
is known. However, it should be noted that the validity of 
this description in the later time of the dynamics is still non- 
trivial. In the noise production simulations from turbulent 
flows, the Lighthill's analogy seems to work well 
Since the noise propagates far from the source region in the 
fluid immediately, complex effects with delayed or nonlin- 
ear interactions seem to be small. 



2.2. Nonlinear feedback in instruments 

In the musical instruments, oscillating energy of the gen- 
erated sound can feed back to the flows after an interaction 
with the resonator. This feedback effect plays the important 
role in the sound generation in musical instruments while it 
can make the problem complicated and difficult ||2l|4||9l. 

In the usual approach to the dynamics in the instruments, 
empirical models are used to simulate the complicated na- 
ture from resonance and oscillation ||4l . In order to clarify 
the mechanism of the oscillation, such a model approach 
often produces good results. However, if we study the de- 
tailed theoretical mechanisms in feedback interactions from 
a sound pressure to the fluids, simulations with more realis- 
tic models should be used. In the next, one of such approach 
using simulation techniques will be presented. 

3. Multiphysics simulation 

Since the acoustic sound is an elastic wave on the air, 
the sound generation by turbulent airflows in musical instru- 
ments is fully simulated by the compressible fluid dynam- 
ics calculation on detailed models of the instruments with 
appropriate boundary and initial conditions. If we require 
sufficient accuracy to the result through the DNS calcula- 
tion, high-performance computing will be necessary. The 
approach of the present work is opposite, where active par- 
titioning of the simulation system into several components 
with different physical properties is introduced. This parti- 
tioning of the whole system is based on the physical insight 
and the integration is performed by the MPS technique. 

In the present problem, the sound generation in musical 
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Figure 2. Mesh structures of the inlet, outlet, 
and resonator. 



instruments, the total system is divided into a part subject 
to incompressible fluid dynamics by Navier-Stokes equa- 
tions and another part for sound propagation and radiation 
by wave equations. These are calculated simultaneously 
with appropriate coupling dynamics between fluid dynam- 
ics and acoustic waves. The schematic picture of this total 
simulation system is given in Fig.[T] 

The model of the instrument in the present work has the 
simplest box-resonator shape (Fig.O since the main subject 
of this work is to demonstrate how useful the multiphysics 
approaches are. The number of generated mesh is shown 
in Table [T] The boundary conditions of this simulation are 
following: At the left inlet, a uniform flow with 1 (m/s) is 
added, and it runs through over the resonator to the right 
outlet where the pressure is given 1 bar (= 10^ Pa). The 
upper wall is a non-friction wall with a velocity 1 (m/s) as 
the fluid, and at the other walls, the velocity is a constant 0, 
and the pressure is zero-gradient. The mesh creation is done 
by blockMesh utility distributed within the open-source 
software package OpenFOAM-1 .4.0 

3.1. Incompressible Large Eddy Simulation 

Flows in the musical instruments are considered as tur- 
bulent while they have relatively low Mach and Reynolds 
numbers. There are many numerical techniques to obtain 
time-evolution of turbulent flows. Recently, Large-Eddy 
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Figure 3. Courant No. of the whole time evo- 
lution (1 second). 



Simulation (LES) is introduced as a rather new tool in the 
field of aeroacoustics, and is often used to simulate compli- 
cated fluid evolutions with turbulence. In the present work, 
LES is used to calculate incompressible flows in the instru- 
ments with appropriate thermo-physical parameters. Al- 
though smaller scales are represented by the so-called sub- 
grid models in this method, it has been studied that acoustic 
sounds generated by turbulent flows can be reproduced by 
the use of this model 1 3 , 8 1 . The LES equations using index 
notation and the summation convention for spatially aver- 
aged field values are given by 



duj ^ dujUj _ _dp_ ^ - ^^'^ 

dt dxj dxi dxj dxj ' 

dxi 

where uj is a velocity field, po is a constant density. 



1 / duj dui 

2 \dxi dxj 



is a strain tensor, and 



T ij — U'iUj U'iUj 



(3) 
(4) 

(5) 

(6) 



^OpenFOAM: http : / . opencf d . co . uk/openf oam/| 



is a subgrid scale tensor. 

The computations are performed on a PC with Dual 
Xeon processors since this is a test execution on the two di- 
mensional minimum model. The numerical scheme is based 
on the finite volume method (FVM) and the time evolution 
scheme is the second order implicit method. The total ex- 
ecution time was about three hours with At = 10~^ (sec). 
The Courant numbers in the whole time evolution are shown 
in Fig. O Snapshots of the results are presented in Fig. IH 
where the absolute values of the velocity u are shown at the 
times t = 0.01, t = 0.1, and t = 1.0 (sec). 




Figure 4. Snapshot of incompressible fluid at t=0.01, 0.1, 1.0 (sec) (from left to right). 



Figure 5. Lighthill's acoustic source term at t=0.01, 0.1, 1.0 (sec) (from left to right). 



3.2. Sound production 

Sound generation and propagation are represented in a 
wave equation with Lighthill's source terms, 



Co 



(7) 



where p' represents a density by sound propagation, and Tij 
is the Lighthill stress tensor which can be approximated by 



Tij = po {UiUj + Tij) . 



(8) 



This term can be obtained even in LES calculation if we use 
the subgrid scale tensor, Eq. ©. 

In Fig. [51 the distribution of the Lighthill's source term, 
the right-hand- side of Eq. (|7]), is shown. From the com- 
parison with Fig. m the sound sources are found around 
small areas where the flows are turbulent and spatially non- 
uniform, i.e., near the corners or edges of walls. In addition 



to these, we can find that vortices become the sound source 
if they move, while it seems that a stable vortex shown in 
the figure of t = 1.0 does not produce sounds. 

3.3. Closing the coupling 

In noise production simulations widely performed by 
DNS and LES, the situation is rather simple since the noise 
is assumed to be spread out. However, in musical instru- 
ments, interactions between airflows and filtered sounds by 
a resonator must be considered. The sound is a density wave 
on the air as a compressible fluid while the MPS system 
studied in the present work is based on the incompressible 
fluid dynamics calculation. In order to take into account 
the feedback effect from the sound to the fluid, extra-terms 
must be introduced. There are no standard and established 
strategies on the feedback interaction so far. 

We can present several forms to introduce the backward 



interaction. If we assume that the air is ideal gas, we can 
use the equation of state for an ideal gas, 



p = pRT 



(9) 



where R is the gas constant and T is the temperature. In 
the incompressible LES calculation, the density p was fixed 
p = Pq. When we define small fluctuation term Ap = 
p' — po by the sound propagated, an additional term 



p^p^ RTAp 



(10) 



can be introduced. Then, this makes a closed loop to the ve- 
locity field through the LES equation Eq. (|3]). What terms 
are necessary and sufficient should be evaluated through 
these simulation studies. 

4. Computational costs 

In this section, we estimate computational costs of mul- 
tiphysics simulations of musical instruments which is com- 
pared to those of the direct numerical simulation. So far in 
the present paper, only two dimensional models are studied. 
In this section, the cost will be estimated for three dimen- 
sional models. 

4.1. Memory requirements 

In our simulations described in the previous section, LES 
calculation scheme is used where the stable results can be 
obtained even if small scaled dynamics arises. However, if 
we use DNS calculation for incompressible fluids based on 
other numerical schemes, finer mesh description is required 
for longer execution of the dynamics. For two dimensional 
models, memory resources required may be still small com- 
pared to the total amount of memory available in the cur- 
rent computers. For three dimensional models, the mem- 
ory increase by representing the fine structures of the com- 
pressible fluids is rather critical and the amount of mem- 
ory required will be an order of lOGBytes. This may not 
be impossible when we implement the parallel computation 
over distributed memory machines, but the performance of 
the program depends deeply on the memory management. 
Thus, the increase of the memory makes longer the execu- 
tion time of the total simulation. 

4.2. Execution time 

When we use the MPS approach, the mesh structures for 
the fluid dynamics and the sound propagation can be defined 
separately. The spatial mesh size Ax for the fluid is deter- 
mined by the finest structure to be studied, and At :^ Ax/U 
of the fluid dynamics can be chosen rather large since the 
velocity U of the fluid is small (~ 1 (m/s) in instruments). 
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Figure 6. Multiscale and multiphysics simula- 
tion describing the nature. 



For representation of the sound propagation, rather coarse 
mesh can be chosen since the spatial structure is not so com- 
plicated around the basic frequency of the resonance. On 
the other hand, time variable must be chosen rather small 
because of the large velocity of the sound (~ 340m/s at nor- 
mal temperature). 

If we perform DNS calculation for the same system, a 
sufficiently fine representation for both the space and time 
must be chosen so that the compressible fluid dynamics with 
the fast density waves is properly simulated. Then, the com- 
putational time will be more than 10^ times longer com- 
pared to the MPS approach. In our two dimensional case, 
the CPU time for one step calculation At = 10~^ is about 1 
second. It seems difficult to execute three dimensional DNS 
calculations unless fine-tuned parallel programs are devel- 
oped and large-scale computer environments are available. 

5. Conclusions 

In the present work, acoustic sound production was sim- 
ulated using MPS approach which couples incompressible 
fluid dynamics calculation by LES method and sound prop- 
agation by wave equations. It was apparent that the amount 
of computational costs is rather small for the MPS case 
compared to the DNS approach, where the separation of 
those programs was essential for systems with different en- 
ergy scales more than 10^. 

This types of multiscale properties can be found in vari- 
ous realistic systems (see Fig. [6]). In the present system, in- 
compressible Navier-Stokes equations are established as ba- 
sic equations to be solved, and the DNS approach still works 
if we could provide sufficient computational resources. On 
the other hand, there are many cases that it is difficult to in- 



troduce a theory of description or a set of equations for the 
whole system, when each component in the system has a 
different scale or the total system is constructed from com- 
ponents with originally different description ifTOl . Descrip- 
tion of the physics is hierarchical, and a special approxima- 
tion valid only to the scale is often used. Then, it is diffi- 
cult to determine which theory (or approximation) is basic. 
Thus, multiphysics or multiscale approaches are significant, 
and direct simulations can be useless even if sufficient com- 
putational resources are available. 
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